Magnetoelastic effects in multiferroic YMnOs 
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We have investigated magnetoelastic effects in multiferroic YMnOs below the antiferromagnetic 
phase transition, Tjv ~ 70 K, using neutron powder diffraction. The a lattice parameter of the 
hexagonal unit cell of YMnOs decreases normally above Tjv, but decreases anomalously below 
Tjv, whereas the c lattice parameter increases with decreasing temperature and then increases 
anomalously below Tjv. The unit cell volume also undergoes an anomalous contraction below 
Tjv. By fitting the background thermal expansion for a non-magnetic lattice with the Einstein- 
Griineisen equation, we determined the lattice strains Aa, Ac and AV due to the magnetoelastic 
effects as a function of temperature. We have also determined the temperature variation of the 
ordered magnetic moment of the Mn ion by fitting the measured Bragg intensities of the nuclear 
and magnetic reflections with the known crystal and magnetic structure models and have established 
that the lattice strain due to the magnetoelastic effect in YMnOa couples with the square of the 
ordered magnetic moment or the square of the order parameter of the antiferromagnetic phase 
transition. 

PACS numbers: 61.05.fm, 65.40.De 



I. INTRODUCTION 



The coupling between spin and lattice degrees of free- 
dom is one of the interesting topics of condensed matter 
researchPJEJ] ft na s drawn renewed interest in connection 
with potentially useful electronic materials, such as colos- 
sal magnetoresistive and multiferroic materials. One of 
the most dominant effects of this coupling is the spon- 
taneous exchange striction associated with the magnetic 
ordering at lower temperatures. If the material does not 
contain any magnetic ions, the temperature variation of 
the volume is smooth and shows no anomalies at low 
temperatures. However, if the material contains mag- 
netic ions that order at low temperature, then at this 
transition temperature the positions of the atoms may be 
modified. The most easily probed effects arc the mod- 
ifications of the unit cell parameters and the unit cell 
volume. Sometimes the lattice symmetry may get modi- 
fied following the symmetry of the magnetically ordered 
state. The symmetry change at the ordering tempera- 
tures also causes crystallographic and/or magnetic do- 
mains to be formed in single crystal samples. In addition 
to the lattice strain effects, important modifications in 
the positional parameters of the atoms and also in the 
bond distances may take place. However, these modifi- 
cations are more difficult to detect because the effects are 
small. To effectively measure this, one needs the temper- 
ature variation of the diffraction intensities and, prefer- 
ably, the diffuse scattering up to high momentum transfer 
of the scattered radiation. A particularly complex situa- 
tion arises for magnetic structures with the propagation 
vector k = 0. In this case, the magnetic reflections occur 
at the same position as the nuclear reflections and the re- 
finement of the magnetic and crystal structures become 
fraught with great difficulties. 



YMn03 belongs to the family of hexagonal manganites 
RMn0 3 (R = Sc, Y, Er, Ho, Tm, Yb, Lu) that show mul- 
tiferroic behaviour^. These hexagonal manganites are 
paraelectric at high temperatures with the centrosym- 
metric space group P&^/mmc. Below about 1000 K 
they undergo a paraelectric-to-ferroelectric transition to 
the non-centrosymmetric structure with the space group 
P63CIT1. Fig. [T] shows schematically the crystal structure 
of YMn(93 in the ferroelectric phase. At lower tempera- 
tures, of the order of about 100 K, the magnetic hexago- 
nal manganites order with a non-collinear antiferromag- 
netic structure with a propagation vector k = 0. Among 
these hexagonal manganites YMJ1O3 and HoMnC>3 have 
been investigated quite intensively^^. Here we report 
the powder neutron diffraction investigation of the tem- 
perature dependence of the crystal and magnetic struc- 
ture of YMnC>3 . We have observed strong magnetoelastic 
effects or spontaneous magnetostriction below the Neel 
temperature Tjv ~ 70 K. The lattice strains are found 
to be proportional to the square of the ordered magnetic 
moments of the Mn ions. 



In order to determine quantitatively the magnetoelas- 
tic effect or the spontaneous magnetostriction, we need to 
measure the temperature variation of the lattice parame- 
ters and the unit cell volume. The high temperature data 
above the ordering temperature can be fitted to either 
a Debye model, using the Griineisen approximation,^ 
or the simpler Einstein model and then extrapolated to 
lower temperatures to give the background variations 
in the absence of magnetism. The simpler Einstein- 
Gruneisen equation serves as an excellent fit function for 
not only the volume but also for the individual lattice 
parameters^ and we have used it here. 
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FIG. 1: (Color online) Schematic representation of the crys- 
tal structure of YMnOa. The red polyhedra represent MnOs 
bipyramids and the blue spheres represent the Y ions. 
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FIG. 2: (Color online) (Upper panel) Temperature variation 
of the diffraction intensities of YMnOs. (Lower panel) The 
same plot in a limited low 26 range showing the temperature 
variation of the peaks with dominant magnetic contribution. 



II. EXPERIMENTAL METHODS 

Neutron diffraction experiments were done on YMn03 
on the high intensity powder diffractometei'2ilD20 of the 
Institute Laue-Langevin in Grenoble. The 115 reflection 
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FIG. 3: (Color online) Rietveld refinement of the (a) crystal 
and magnetic structures of YMnOs at T = 2 K and (b) the 
crystal structure of YMn0 3 at T = 300 K. The ticks in (a) 
give the positions of the nuclear and magnetic Bragg peaks 
whereas those in (b) give the positions of the nuclear peaks 
only. 



from a Ge monochromator at a high take-off angle of 118° 
gave a neutron wavelength of 1.868 A. Approximately 3 g 
of powdered YMn03 was placed inside an 8 mm diameter 
vanadium can, which was fixed to the sample stick of a 
standard 4 He cryostat. 



III. RESULTS 

We measured the diffraction intensities from YMnOs 
as a temperature in the range 2 — 300 K. The upper- 
panel of Fig. [2] shows the temperature variation of the 
diffraction intensities of YMn03 whereas the lower panel 
shows a similar plot in a limited low 29 range showing 
the temperature variation of the peaks with a dominant 
magnetic contribution. The magnetic contributions to 
the intensity are clearly seen in the plot of the lower 
panel. The reflections like hOl with 1 = odd disappear 
completely at T/v as they are forbidden reflections in the 
space group PG^cm and have only magnetic contribution 
and therefore become zero at T/v. The other reflections 
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FIG. 4: (Color online) Temperature variation of the lattice 
parameters a, c, and the unit cell volume V of YMnOs plotted 
on the left panel. The red curves in these figures represent 
the lattice parameter and the unit cell volume obtained by 
fitting the high temperature data by the Einstein-Griineisen 
equation and extrapolated to low temperature to give the 
background for the non-magnetic lattice. On the right panel 
the temperature dependence of the lattice strains Aa, Ac and 
AV obtained by subtracting the non-magnetic background 
have been plotted. 



that have both magnetic and nuclear contributions at- 
tain a constant value (nuclear) above Tjv. The Rietveld 
refinement!^ of the diffraction data was done using the 
Fullprof prograrrPHl. The refinement results from YMJ1O3 
at T = 2 K and 300 K are shown in Fig. |) At T = 2 K, 
the intensity data have been fitted with the crystal and 
the magnetic structure modeP2l f YMnC>3. At T = 300 
K, in the paramagnetic state, the refinement was done 
with the crystal structure model only. 

Figure [4] shows the temperature variation of the lattice 
parameters a, c, and the unit cell volume V of YMJ1O3 
plotted on the left panel. The lattice parameter a de- 
creases continuously with decreasing temperature but, 
close to the Neel temperature Tjv ~ 70 K, shows a mag- 
netoelastic or magnetostriction anomaly. The red curve 
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FIG. 5: (Color online) Temperature variation of the ordered 
x-component of magnetic moment m x of the Mn ion obtained 
by fitting the measured Bragg intensities of the nuclear and 
magnetic reflections with the known crystal and magnetic 
structure models. 



represents the background variation of the a lattice pa- 
rameter for a non-magnetic solid obtained by fitting the 
data in the paramagnetic state by the Einstein-Griineisen 
equation. By subtracting the background from the data 
we determined the lattice strain Aa, plotted on the right 
panel of Figure [4] The c lattice parameter increases 
continuously with decreasing temperature down to 7jv 
and then shows an anomalous increase or positive mag- 
netostriction. By subtracting the background variation 
of the c lattice parameter (red continuous curve) from 
the data we determined the lattice strain Ac plotted 
on the right panel. Similar plots for the unit cell vol- 
ume V and and the volume strain AV are shown in 
Figure [4] in the left and the right panels, respectively. 
The c/a ratio, which is 1.8572 ± 0.0001 at T = 296 K, 
increases linearly with decreasing temperature down to 
the Neel temperature and then attains a constant value 
c/a = 1.8643 ± 0.0001 at low temperatures. A closer 
inspection of the temperature variation of the lattice pa- 
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FIG. 6: (Color online) Linear plots of the lattice strains Aa, 
Ac and AV vs. the square of the ordered magnetic moment 
of the Mn ion in YMn0 3 . 



rameters shows that the magnetoelastic effects become 
actually noticeable at temperatures substantially above 
TV rj 70 K. We can tentatively assign a temperature 
T* « 100 K at which the magnetoelastic effect first be- 
comes noticeable. The appearance of the magnetoelas- 
tic effect at temperature higher than Tjy is presumably 
due to the short-range spin correlations. The sponta- 
neous linear magnetostriction along the a-axis at T = 
is Aa/ao = —2.944 x 10~ 4 and that along the c- 
axis is Ac/cq — 2.633 x 10~ 4 . The spontaneous vol- 
ume magnetostriction AV/Vq is —3.252 x 10~ 4 . Here 
ao, Co and Vo are the extrapolated lattice parameters 
and the unit cell volumes at T = for a fictitious non- 
magnetic solid obtained by the method described in the 
previous section. Note that although we determined 
AV/Vq = -3.252 x 10~ 4 by fitting the non-magnetic 
background and subtracting the experimentally deter- 
mined unit cell volume, one can also calculate this from 
the relation AV/Vq = 2 Aa/ao + Ac/co- The calcu- 
lated value turns out to be AV/Vq — —3.255 x 10~ 4 , 
very close to the value —3.252 x 10 -4 determined semi- 



independently. This indicates that the method we used 
for the determination of the non-magnetic background is 
reliable or at least consistent. One must however point 
out that the paramagnetic and fictitious non-magnetic 
background volumes are not really the same but we have 
assumed this to be the case for simplicity during the 
present analysis. 

Figure [5] shows the temperature variation of the x- 
component of the ordered magnetic moment m x of the 
Mn ion, obtained by fitting the measured Bragg inten- 
sities of the nuclear and magnetic reflections with the 
known crystal and magnetic structure models. The or- 
dered magnetic moment decreases with increasing tem- 
perature and becomes zero at Tjy ~ 75 K. We did not 
have enough data points close to T/v and therefore the 
critical exponent (3 could not be determined. 

One expects that the spontaneous magnetostriction is 
proportional to the square of the ordered magnetic mo- 
ment. This is at least the case for volume magnetostric- 
tion for a single sublattice ferromagnelP^. However our 
recent investigations^^!] on transition metal difluorides 
MF 2 (M = Mn, Co, Fe, Ni) showed that this is also the 
case for these simple antiferromagnets with the rutile 
structure. To check the coupling between the sponta- 
neous magnetostriction and the ordered magnetic mo- 
ment, or the order parameter of the antiferromagnetic 
phase transition in YMn03, we plotted in Fig. [6] the 
lattice strains Aa, Ac and AV vs. the square of the x- 
component of the ordered magnetic moment of the Mn 
ion. The linear plots show that the lattice strains in 
YMn03 indeed couple to the square of the magnetic mo- 
ment or the order parameter. This result also suggests 
that a second-order spin Hamiltonian is sufficient to de- 
scribe the magnetic system of YMn03. It is not necessary 
to invoke higher order spin Hamiltonian consisting of, for 
example, four-spin exchange for YM11O3. 



IV. DISCUSSION AND CONCLUSIONS 

Although we were successful in determining the mag- 
netoelastic effects on the lattice parameters, our attempts 
to extract the temperature variation of the positional pa- 
rameters of YMn03 below T/v did not succeed, presum- 
ably because of the k = magnetic structure of YMn03. 
In such structures, magnetic reflections arise below TV on 
the same positions as the reflections of structural origin. 
Well above T/v diffuse magnetic scattering is observed, 
becoming magnetic Bragg reflections below T/y. How- 
ever, the magnetic scattering in general coincide with 
nuclear reflections. There are only a few pure magnetic 
reflections where structural reflections are not allowed 
due the space group extinctions. The least squares re- 
finement of the crystal and magnetic structures using 
overlapping nuclear and magnetic intensities at the same 
reciprocal lattice positions has the problem of separat- 
ing magnetic and nuclear intensities and therefore lead 
to correlations between the magnetic and structural pa- 
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rameters. Due to the correlations between structural and 
magnetic parameters, the positional parameters show ab- 
normal deviations at and near Tjy. Also these deviations 
are different if you come from the high or low tempera- 
ture sides in a sequential refinement in which one uses the 
refined parameters of the previous data set as the start 
parameters of that for the next data set. One way to 
avoid such unwanted correlations is to use polarized neu- 
tron diffraction. The other useful method to disentangle 
the magnetic and structural parameters is to measure 
neutron diffraction intensities to a very high momentum 
transfer Q = iir^j^-. Our present data collected on a re- 
actor based diffractometer at the thermal neutron beam 
are limited to about 7 A~ x only. Such high-Q diffrac- 
tion data can be collected at a spallation neutron source. 
One can then refine the crystal structure using only the 
high-Q data and the magnetic structure using only the 
low-Q data. In this way one can hope to get the magne- 
toelastic effect on the positional parameters (bond dis- 
tances, symmetric shifts etc.) as well and not just the 
lattice magnetoelastic effect that we determined in the 
present paper. Lee et alP^have obtained magnetoelastic 
effects of the structural parameters and bond distances 
of YMJ1O3, LuMn03 and Y0.5LU0.5MJ1O3 from their neu- 
tron and X-ray diffraction investigations but they have 
not explained how they have avoided such correlations 
between the structural and magnetic parameters. In the 
absence of such explanation it is difficult to accept the 
results given by these authors^. It is extremely impor- 
tant to check the correlation matrix of the refined struc- 
tural and magnetic parameters, especially in the case of 



a k = magnetic structure, because the structural and 
magnetic intensities appear at the same reciprocal lattice 
points and least squares refinement process involves di- 
viding the intensities into structural and magnetic parts 
rather arbitrarily. Only polarized neutron diffraction can 
separate the magnetic and structural contributions di- 
rectly and the authors^ did not use such a method. We 
also checked carefully the supplementary information of 
the paper in the on-line version. It contains a lot of other 
information but does not contain the crucial information 
about the correlation matrix of the least squares refine- 
ment. 

In conclusion, we have investigated the spontaneous 
magnetoelastic effects or magnetostrction in multiferroic 
YM11O3 by neutron powder diffraction. By subtract- 
ing the background thermal expansion for a nonmag- 
netic lattice from the measured thermal expansions we 
determined the lattice strains Aa, Ac and AV^ due to 
the magnetoelastic effects as a function of temperature. 
We have also determined the temperature variation of 
the ordered magnetic moment of Mn ion by fitting the 
measured Bragg intensities of the nuclear and magnetic 
reflections with the known crystal and magnetic struc- 
ture models. Our results show that the lattice strain due 
to the spontaneous magnetostrictions in YMn03 couples 
with the square of the ordered magnetic moment or the 
square of the order parameter of the antiferromagnetic 
phase transition. This result also suggests that a second- 
order spin Hamiltonian is sufficient to describe the mag- 
netic system of YMn03. 
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